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SUMMARY 


We  have  begun  the  calculation  of  relaxation  about  the  native  point  defect  in  HgCdTe. 
We  have  found  that  approximate  self-consistent  schemes  are  inadequate  for  the  large 
displacements  of  the  highly  strained  defects;  therefore,  these  calculations  are  being  done  self- 
consistently,  and  thus  are  computationally  time-consuming.  The  fields  produced  by  long-range 
strain  fields  of  dislocations  via  the  piezoelectric  effect  and  charged  dislocation  cores  have  been 
calculated.  We  find  that  the  piezoelectric  potential  is  a  superposition  of  cylindrical  quadrupole 
and  hexadecapole  terms.  The  impact  of  dislocation  core  charges  depends  on  the  location  of  the 
dislocation,  and  is  greater  in  the  depletion  region  where  the  screening  is  reduced,  and  smaller  in 
the  neutral  regions  of  the  device.  We  also  find  that  the  impact  decreases  as  the  temperature  is 
increased.  We  have  completed  the  calculations  of  the  unrelaxed  native  point  defect  energies  in 
ZnSe  within  the  local  density  approximation  and  have  begun  the  calculation  of  the  gradient 
correction  to  these  energies.  Based  on  these  preliminary  numbers,  we  predict  the  zinc  vacancy 
and  the  zinc  antisite,  and  the  zinc  interstitial  to  be  the  dominant  defects  in  ZnSe.  We  are 
currently  calculating  the  relaxation  about  the  most  important  defects  in  ZnSe. 
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1. 


NATIVE  POINT  DEFECTS  IN  IRFPA  ACTIVE  AND  SUBSTRATE  MATERIALS 


1 .1  NATIVE  POINT  DEFECTS 

In  the  first  quarter,  we  reported  the  results  of  native  defect  concentrations  in 
Hgo.8Cdo.2Te.  During  the  second  quarter,  we  have  focused  on  refining  the  defect  formation 
energies  used  in  those  predictions.  In  particular,  we  have  begun  the  calculation  of  the  relaxation 
energies  about  the  native  defects  in  HgCdTe.  We  are  focusing  on  the  important  defects  in 
Hgo.8Cdo.2Te,  and  the  mercury  interstitial,  which  is  important  to  mercury  diffusion  and  for 
which  strain  energies  are  expected  to  be  large. 

Our  original  scheme  to  calculate  the  relaxations  using  the  Harris  approximation,  in  which 
the  position  of  the  relaxed  atoms  is  calculated  non-self-consistently,  has  proved  to  be  inadequate 
for  cases  in  which  the  relaxations  are  large.  Thus,  the  position  of  the  relaxed  atoms  must  be 
found  self-consistently.  We  are  currendy  completing  these  self-consistent  calculations  for  the 
mercury  vacancy,  the  tellurium  antisite,  and  the  mercury  interstitial.  We  expect  these 
calculations  to  be  completed  in  the  next  quarter  and  improved  predictions  of  defect 
concentrations  to  be  made  at  that  point.  The  relaxation  calculations  are  also  being  extended  to 
larger  supercells.  These  calculations  are  computationally  time-consuming,  and  availability  of 
supercomputer  time  would  expedite  our  work. 

1 .2  DISLOCATION  STRAIN  FIELDS  AND  PIEZOELECTRIC  POTENTIALS 

We  pointed  out  in  Quarterly  Technical  Report  1  that  the  MCT  crystal  structure  lacks 
inversion  symmetry  and  therefore  is  piezoelectric,  and  that  the  strain-induced  bound  charge 
density  and  the  fields  that  it  generates  can  affect  the  transport  properties  of  the  semiconductor. 
We  reported  on  a  first  approach  to  this  issue,  and  we  have  continued  the  work  in  the  current 
quarter. 

The  strain  field  around  a  dislocation  can  generate  an  electric  field,  and  conversely  an 
electric  field  applied  to  a  piezoelectric  crystal  can  produce  a  stress  that  may  deform  the  crystal. 
This  self-consistent  field  problem  has  been  treated  rigorously  for  a  straight  dislocation  in  a 
homogeneous  insulating  crystal  (intrinsic  wide-gap  semiconductor)  by  Saada  (1971).  The 
general  analysis  was  extended  by  Faivre  and  Saada  (1972)  to  include  the  effects  of  screening  by 
free  carriers.  Several  qualitatively  significant  conclusions  can  be  drawn  from  this  work:  ( 1 )  The 
strains  are  only  slightly  affected  by  the  electric  field;  that  is  to  say,  the  inverse  piezoelectric 
effect  is  small,  and  the  electric  field  can  be  calculated  directly  without  imposing  the  requirement 
of  self-consistency.  (2)  The  electric  field  in  the  vicinity  of  the  dislocation  core  can  be  quite 
large,  and  the  dislocation  can  interact  strongly  with  charges  (point  defects  or  mobile  carriers)  in 
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its  vicinity.  (3)  In  a  semiconductor  the  dislocation  preferentially  traps  carriers  of  a  panicular 
sign.  (4)  The  estimated  screening  distance  is  of  the  order  of  several  Debye  lengths. 

Booyens  and  Vermaak  (1979)  have  applied  Saada’s  analysis  to  (1 10)  dislocations  in 
zincblende  III-V  semiconductors,  and  Booyens  and  Basson  (1984)  to  CdTe  and  MCT.  They  do 
not  explicitly  include  the  effects  of  screening  as  analyzed  by  Faivre  and  Saada  (1972)  in  these 
calculations.  As  a  result,  they  find  that  an  uncharged  dislocation,  in  the  absence  of  screening, 
generates  a  net  radial  field.  They  suggest  that  the  dislocation  core  is  likely  to  acquire  a  charge 
that  compensates  for  the  radial  field  and  so  minimizes  the  total  energy  of  the  dislocation.  While 
this  is  a  plausible  suggestion,  it  does  depend  on  unknown  details  of  the  core  electronic  states  and 
needs  to  be  examined  more  closely.  We  have  similarly  extended  the  analysis  to  the  commonly 
observed  60°  dislocations  in  intrinsic  MCT  and  obtained  detailed  plots  of  the  electric  potential 
distribution  arising  from  the  strain  about  a  straight  dislocation.  We  find  that  the  potential  is  a 
superposition  of  cylindrical  quadrupole  and  hexadecapole  terms.  The  theory  and  initial  results  of 
the  calculation  are  given  in  the  Appendix.  These  results  are  preliminary,  since  they  do  not  yet 
take  into  account  screening  by  free  carriers,  anisotropy,  and  interactions  among  a  distribution  of 
dislocations. 

1 .3  DISLOCATION  CORE  CHARGES 

The  fields  associated  with  electrically  charged  dislocation  cores  can  act  as  scattering  and 
recombination  centers  and  may  attract  or  repel  charged  point  defects.  As  pointed  out  above  in 
connection  with  the  strain-induced  charges  and  fields,  the  presence  and  amount  of  charge  on  a 
dislocation  core  depend  on  the  distribution  of  gap  levels  in  the  core  and  on  the  Fermi  level  in  the 
semiconductor.  Here  we  take  a  phenomenological  view  and  consider  the  electric  field  and 
potential  generated  by  a  given  core  charge,  determined  from  the  band  bending  around  an  inverted 
core.  The  band  gap  of  the  material  sets  the  boundary  condition  for  the  potential  at  the  core 
radius.  We  assume,  conservatively,  a  core  radius  of  5  A  and  a  core  edge  potential  of  0. 1  eV. 
equal  to  the  band  gap.  The  effective  radius  of  the  field  is  taken  to  be  that  distance  from  the 
dislocation  core  within  which  the  field  is  large  enough  to  move  a  carrier  the  distance  between 
dislocations  in  a  minority  carrier  lifetime  (as  measured  in  dislocation-free  material).  The 
rationale  for  this  assumption  is  that  a  field  smaller  than  this  is  not  likely  to  affect  the  carrier 
lifetime  or  mobility.  If  we  use  a  lifetime  of  10’^  s  and  a  mobility  of  10^  cm^/Vs,  this  sets  the 
minimum  effective  field  at  approximately  lOno'*^  V/cm,  with  the  dislocation  density  n[)  in  cm  -. 
This  is  used  to  determine  a  numerical  threshold  field. 

The  principle  of  the  calculation,  both  in  the  depletion  approximation  and  with  screening, 
was  given  in  Appendix  B  of  Quarterly  Report  1 .  We  have  now  extended  the  work,  using  a  better 
value  of  the  optical  range  dielectric  constant  (Brice  and  Capper  1987),  and  taking  into  account 
the  temperature  dependence  of  the  Debye  frequency. 

The  results  are  summarized  in  Figure  1.  The  figure  shows  two  sets  of  results:  the  radius 
R  of  the  field  region  surrounding  the  positively  charged  core  as  a  function  of  carrier  density 
calculated  on  the  basis  of  the  depletion  approximation,  and  the  screening  radius  /  of  the  region 
inside  which  the  electric  field  due  to  the  charged  core  exceeds  10  mV.  This  threshold  for  the 
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10lO  10^^  10''2  10^^  10^'*  10^5 


carrier  density  n,  cm^ 

Figure  1 .  Dependence  of  space  charge  radius  R  (depletion  approximation) 
and  screening  radius  t  on  carrier  density.  For  the  range  shown 
here,  we  find?  a 

be-12089:i-sd 

electric  field  is  based  on  the  experimental  observations  of  Johnson  et  al.  (1992)  and  Shin  et  al. 
(1992)  that  the  RqA  values  and  minority  carrier  lifetimes  decrease  rapidly  with  the  dislocation 
density  no  when  it  exceeds  10^  cm'2.  TTie  calculation  of  the  screening  radius  i  is  more  reliable 
than  the  depletion  approximation,  since  it  takes  into  account  the  thermal  agitation  of  the  earners. 
With  the  parameters  we  use,  the  calculation  of  the  screening  radius  over  five  orders  of  magnitude 
of  the  carrier  density,  and  over  the  temperature  range  of  40  to  300  K,  results  in  values  of  f  that 
range  from  about  5  Lp  (Debye  lengths)  at  n  =  10*®  cm'^  to  10  Lp  at  n  =  10*^  cm’^,  resulting  in 
an  overall  decrease  of  ^  with  carrier  density  at  a  rate  slightly  slower  than 

The  carrier  density  in  the  neutral  region  of  the  extrinsic  device  material  generally  equals 
or  exceeds  10*^  cm*^.  The  screening  length  at  this  carrier  density  even  at  room  temperature  is 
less  than  4  pm,  and  only  about  1.5  pm  at  40  K.  Since  the  mobile  carrier  density  in  the  neutral 
regions  generally  exceeds  10*"*  cm*^,  we  can  conclude  that  in  the  presence  of  mobile  earners, 
especially  at  cryogenic  temperatures,  the  core  charge  fields  have  a  limited  range,  in  rough 
agreement  with  the  observed  small  variation  of  the  responsivity  with  dislocation  density. 


The  situation  is  less  clear  in  an  extended  depletion  region,  such  as  the  vicinity  of  the 
metallurgical  junction  of  the  detector.  In  this  region  the  carrier  density  is  reduced  to  the  intrinsic 
value  nj  under  zero  bias,  and  below  this  with  a  reverse  bias.  We  estimate  a  value  of  ni  of  about 
5x  cm‘3  for  Hgo.8Cdo.2Te  at  40  K,  with  a  decrease  of  approximately  a  factor  of  2  for  every 
5  mV  of  reverse  bias  caused  by  the  separation  of  the  quasi  Fermi  levels.  This  reduction  of  the 
carrier  density  is  confined  to  a  small  region,  and  an  assessment  of  its  quantitative  effect  will 
require  detailed  computation.  Nonetheless,  it  is  worth  observing  that  for  this  carrier  density  we 
estimate  a  screening  radius  of  17  pm,  consistent  with  the  experimental  finding  (Johnson  et  al., 
1992)  that  RoA  at  this  temperature  is  already  decreasing  with  dislocation  density  in  this  range. 
The  literature  values  of  nj  =  1 .2  x  10'^  cm-3  at  77  K  and  ni  =  4.0  x  10*^  cm-^  at  300  K  (quoted  in 
Brice  and  Capper,  1987)  correspond  to  progressively  shorter  screening  lengths,  which  is  also 
consistent  with  the  experimental  observation  that  the  decrease  in  RoA  with  dislocation  density 
has  a  higher  threshold  at  the  higher  temperature. 

Thus,  to  summarize  we  find  that  dislocations  impact  device  performance  through  two 
possible  mechanisms:  core  line  charge  accompanied  by  a  depletion  layer  with  radial  field 
surrounding  it,  and  an  angular  potential  variation  arising  from  the  long-range  strain  and  the 
piezoelectric  effect.  In  a  neutral  region  of  the  semiconductor — that  is,  the  n-region  below  the 
depletion  layer  in  a  p-on-n  heterojunction  diode — screening  prevents  the  dislocations  from 
greatly  influencing  device  characteristics.  As  a  consequence,  the  responsivity  of  such  devices, 
which  is  dominated  by  the  minority  carrier  recombination  lifetime,  varies  slowly  with  dislocation 
density.  However,  in  the  depletion  region  screening  is  greatly  reduced,  and  the  long-range 
angular  potential  produced  by  a  dislocation  that  threads  through  the  depletion  region  is  capable, 
in  principle,  of  increasing  the  tunneling  currents  at  low  temperature  as  observed.  We  need  an 
accurate  value  for  the  HgCdTe  piezoelectric  coefficient  and  a  determination  of  the  core  charge 
states  of  the  various  dislocation  types  before  we  can,  with  more  assurance,  specify  the 
mechanism  through  which  dislocations  affect  the  responsivity  and  noise  of  devices.  Some 
dislocation  types  are  certain  to  be  more  troublesome  than  others  in  a  given  device  structure. 

Once  the  worst  offenders  are  identified  we  can  begin  to  design  ways  to  eliminate  or  passivate 
them. 
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2.  PHOTONIC  MATERIALS 


2.1  WIDE-GAP  ll-VI  COMPOUNDS  (ZnSe  AS  PROTOTYPE) 

We  have  completed  the  local  density  calculation  of  the  fomiation  energies  for  the  neutral 
unrelaxed  native  defects  in  ZnSe.  The  defects  we  have  considered  are  the  cation  and  anion 
vacancies  Vzn  and  Vse,  the  cation  and  anion  antisite  Sezn  and  Znse  and  the  cation  and  anion 
interstitial  Sej  and  Zni.  The  interstitial  site  considered  is  tetrahedral,  where  the  four  first 
neighbors  are  cations;  we  will  later  return  to  calculation  for  the  tetrahedral  site  where  the  four 
first  neighbors  are  anions,  '^he  defect  reactions  considered  are 


ZnSe  — »  VznSe  -t-  Zng 

(1) 

Zng  ->  Zn  Vse 

(2) 

2Zng  Zn  Znse 

(3) 

2ZnSe  ->  SeznSe  -t-  2Zng 

(4) 

Zng  — >  Znj 

(5) 

ZnSe  — >  Sej  +  ^tig 

(6) 

where  Zng  is  the  zinc  in  the  free  atom  state. 

The  calculation  of  the  defect  formation  energies  consists  of  several  steps.  First,  the 
formation  energy  from  the  unrelaxed  lattice  is  calculated.  These  energies  for  the  16-atom 
supercell  are  given  in  Table  1.  To  these  energies  the  gradient  correction  to  the  local  density 
formation  energies  must  be  added.  This  correction  is  important  when  the  free  atom  is  used  as  a 
reference  state,  as  is  the  case  when  we  wish  to  consider  the  vapor  pressure  as  an  external  free 
variable.  The  gradient  corrections  for  the  defects  have  been  completed  using  a  coarse  mesh  for 
the  k-space  Brillouin  zone  integration.  Based  on  these  calculations,  the  gradient  corrections  for 
defect  reactions  1  through  6  are  sizable,  ranging  from  0.2  to  0.8  eV,  and  thus  merit  a  more 
precise  calculation.  The  self-consistent  charge  densities  from  the  coarse  mesh  calculation  are 
being  used  to  calculate  the  gradient  correction  self-consistently  on  a  finer  mesh.  1  hese 
calculations  are  nearly  done  and  will  be  reported  in  the  next  quarter. 

The  relaxation  energy  must  be  calculated  by  allowing  the  overall  lattice  constant  to  obtain 
a  minimum,  and  searching  for  the  minimum  energy  position  of  the  near-neighbor  atoms.  The 
former  relaxation  has  been  completed,  and  this  relaxation  energy  is  included  in  the  energies 
given  in  Table  1. 
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Table  1 


DEFECT  FORMATION  ENERGIES  AND  ENTROPIES  FOR  THE  NEUTRAL,  UNRELAXED  DEFECTS. 

Energies  do  not  Include  the  gradient  correction 


Defect 

Energy  (eV) 

Entropy  (ke) 

VZn 

5.00 

-7.58 

VSe 

0.76 

8.34 

Znse 

-1.23 

15.9 

Sezn 

9.40 

-16.0 

Zni 

1.08 

7.91 

Sei 

7.37 

-7.91 

We  are  currently  calculating  the  second  contribution  to  the  relaxation  energy.  We  have 
been  experimenting  with  non-self-consistent  schemes  whereby  the  charge  densities  about  each 
atom,  for  example  in  the  unrelaxed  configuration,  are  used  to  calculate  the  energy  for  the  lattice 
in  a  relaxed  configuration  using  a  Harris  approach.  While  this  method  appears  to  be  reliable  for 
small  distortions,  it  appears  to  be  unreliable  for  the  large  distortions  we  are  finding  about  the 
interstitial  and  antisite  defects.  Thus,  we  are  forced  to  complete  self-consistent  calculations  of 
the  relaxation,  which  entail  a  large  amount  of  computer  time.  Because  our  rate  of  progress  is 
limited  by  the  speed  of  our  in-house  computer,  which  is  timeshared  among  several  projects, 
availability  of  supercomputer  resources  from  DARPA  would  expedite  these  calculations.  Our 
codes  have  previously  been  optimized  for  operation  on  many  machines,  including  the  CRAY 
series. 

Temperature-dependent  defect  formation  entropies  have  been  calculated.  These  entropies 
are  a  necessary  ingredient  in  the  calculation  of  the  native  defect  densities.  When  a  defect  is 
created  in  the  ZnSe  lattice,  the  phonon  modes  of  the  crystal  are  modified,  and  these 
modifications  contribute  to  the  defect  formation  entropy.  To  calculate  this  contribution  to  the 
defect  formation  entropy,  we  have  used  a  valence  force-field  model  for  the  lattice  phonons.  A 
Green’s  function  method  is  then  used  to  introduce  a  localized  perturbation,  for  example  a 
vacancy  or  an  antisite,  into  the  lattice.  From  the  change  in  the  density  of  (vibrational)  states,  the 
defect  formation  entropy  can  be  calculated.  The  formation  entropies  in  units  of  Boltzmann’s 
constant  for  800  °C  are  given  in  Table  1 .  These  entropies  correspond  to  the  reactions  in  Eqs.  1 
through  6  above. 

2.2  NONLINEAR  OPTICAL  MATERIALS  (LINbOs  AS  PROTOTYPE) 

The  work  on  LiNbOs  is  scheduled  to  start  in  the  next  quarter. 
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3.  WORK  PLANNED 


During  the  next  quarter,  we  expect  to  complete  the  calculation  of  the  relaxation  of  the 
near  neighbors  about  the  most  important  native  defects  in  both  HgCdTe  and  ZnSe,  and  to 
complete  the  gradient  correction  to  the  defect  energies  in  ZnSe.  We  expect  to  make  preliminary 
estimates  of  native  defect  densities  in  ZnSe,  and  to  refine  those  in  HgCdTe.  Work  will 
commence  on  the  properties  of  LiNb03. 
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Appendix 


PIEZOELECTRIC  POTENTIAL  AT  A  60*  DISLOCATION  IN  SPHALERITE 


The  dislocation  lies  along  [lOl]  and  has  a  Burgers  vector  b  at  60°  to  the  dislocation  line 
The  edge  component  of  the  Burgers  vector  b®  is  chosen  as  xi.  X2  is  the  normal  to  the  slip 
plane,  and  X3  is  ^  (Figure  A-1).  The  transformation  matrix  from  principal  axis  coordinates  is 


f  1-2  1  " 

^  V2>(rv2 


The  transformation  of  the  piezoelectric  tensor  is 

d  ijk  =  TjpTjqTkjdpqr. 

In  cubic  symmetry  and  principal  axis  coordinates  all  the  nonzero  components  of  the 
piezoelectric  tensor  are  equal  di23=di32=d23i=d2i3=d3i2=d32i=d. 


The  nonzero  components  of  the  transformed  tensor  (d’)  are 

2 

di33  =  dsn  =  dssi  =  “dm  “ 


dl  12  -  dl21  -  d21 1  -  d233  -  d323  “  d332  “ 


_2_ 

V3' 


d222  -~<l  ’ 


in  two-index  (Voigt)  notation 


(d')  = 


(11  1  ^ 

— =d  0  -=d  0  0  — =d 

V6  V6  V3 

— 5=d-^d— ^d  0  0  0 

V3  VJ  V3 


V 


0  0  0  — ^d  4rd  0 

^  V6 


The  piezoelectrically  induced  bound  charge  density  is 

a 


♦  - « 


(1) 


1 


(a)  Definition  of  coordinate  system 


/6 


(b)  Orientation  of  line  vector  ^  and  Burgers  vector  b  of  a  60°  dislocation 


Figure  A-1 .  Geometry  used  in  calculating  properties  of  dislocations 


bc-220893-!,d 


where  e  is  the  isotropic  dielectric  constant,  ejk  is  the  jk  component  of  the  strain  tensor,  and 
summation  over  repeated  indexes  is  implied. 


For  a  screw  dislocation  in  linear  anisotropic  theory  the  displacements  are 

.Ui  =  U2  =  0  , 


U3  =  -  —  tan  'a 

2jt 


^2^ 


xi; 


where 


A  =  c'44  =  C44-tH;  C55  =  C44-7H:  H  =  2C44  +  C.2-C..  . 

C44  3  6 

Here,  Cij  and  c'jj  are  the  components  of  the  elastic  tensor  in  the  principal  axis  and 
transformed  coordinate  systems,  respectively.  The  nonzero  components  of  the  strain  tensor  are. 
using 

ejj  =  y  (uj  j  +  Uj  j):  ei3  =  e3)  and  e23  =  e32;  all  others  zero. 

Using  (1)  and  (2)  shows  immediately  that  for  this  case  p  =  0. 

For  an  edge  dislocation,  with  v  Poisson’s  ratio 


-if  ^2 


tan  — 
xi 


XiX2 


2(l-v)(x^+x|)J 


l-2v 
4(1 -v) 


ln(x^+X2)  + 


2 

X1-X2 


4(1-v)(xVx2) 


U3  =  0. 


The  nonzero  strain  components  are 

b  (1-2v)x2+(3-2v)xjX2 

2(v-l)(x^-»-x|)^ 

b  (1-2v)x2-(1+2v)x^X2 

2(v-1)(x^+X2)^ 

.  2  3 

b  X1X2-X1 

2(v-l)(x^+xl)^  ' 
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From  (2)  and  ( 1 ) 


p(r)  — (diiieii.i  +  2dii2ej2,i  +  d2iien.2  +  d222  e22.2) 


so  we  need  the  derivatives 


b  (2v+1)xjX2+(2v-3)x]x2 

b  (2v-l)x2-6xjx|+(3-2v)xj 
2(v-l)(xi+x|)^ 


®22.2  - 

®12.1  = 


b  (2v-l)x2+6xjx|-(2v+l)xi 

2(v-l)(x^i+x|)^ 

4  ,22,  4 
b  X2-6X;X2+Xi 

2(v-l)(xVx^)^ 


giving 


p(xi  .X2)  =  Cl  1.1  -  2ei2,i  -  eii,2  +  2e22.2) 


bd  (2v  -  3)x;-2V2(2v  +  1)x,x|h-30x  iX;+2V2(3  -  2v)x^x,-(2v  +  7)x^ 
2n^J3t  2{v-1)|x^+X2)^ 


When  this  expression  is  converted  to  polar  coordinates  in  the  xix2  plane  (xyplane)  this 
becomes  a  somewhat  simpler  expression  that  has  a  radial  dependence  with  angular  terms 
having  quadrupole  and  hexadecapole  form: 


P  = 


bd 

27C'/3e 


(P1  +  P2+P3  +  P4  +  P5) 
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with 


(2v-3)  y'*  (2v-3)  cos40-4cos20+3 

"  2(v-l) 

_  -2V2(2v+l)  xy^  _  -2V2(2v-H)  2sin2e-sin49 
^2“  2(v-l)  "  2(v-l)  gr2 

_  30  x^y^  _  30  l-cos4e 

‘^3"2(v-1)  “2(v-l)  gr2 

2V2(3-2v)  x\  2V2(3-2v)  sin4e+2sin2e 
2(l-v)  ■  2(v-l)  gr2 

_  -(2v+7)  _  -(2v+7)  cos49+4cos204-3 

2(v-l)  ~  2(v-l)  gr2 


p  = - - -r  [5cos46-V2  sin46+-2(l+v)cos26-V2  (l-v)sin20l 

4ttV3e(l-v) 


=  -y  ^  [EnCosCnO)  +bnSin(n0) 
r  n=2,4 

To  obtain  the  piezoelectric  potential  we  need  to  solve  the  Poisson  equation 

v2<l>  =  -ilP 

e 


with  this  charge  distribution  as  the  source.  In  two  dimensions  the  formal  solution  of  the  Poisson 
equation  is 


In— + 


m 

cosm(0-0') 


p(r') 


d0'[sin(m0)sin(m0’>+cos(m0)cos(m0')lx 


[ancos(n0')+b_sin(n0’)]  f  — ^ 

■^2.4  ■’» 


vr> 


where  r>  and  r<  are  the  greater  and  lesser,  respectively,  of  r  and  r'. 
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The  integral  diverges  at  the  origin  because  of  the  dependence  of  p.  To  deal  with  this 
difficulty  we  adopt  the  same  procedure  that  is  used  to  avoid  the  divergence  of  the  elastic  energy 
of  a  dislocation  at  the  origin  in  continuum  theory.  We  set  cutoffs  at  an  inner  radius  Ri  and  an 
outer  radius  R2.  The  cutoff  at  Ri  removes  the  core  singularity,  and  the  cutoff  at  R2  allows  p  to 
be  expanded  in  multipoles.  Later  we  will  identify  Ri  as  the  core  radius  and  we  will  let  R2  —> 
We  now  have 


<I>(r,0)  =  — anCOs(n0)+bnSin(n6)[  —  f  (r')"~*dr' +  r"f  (r')  "  'df 

Vr"J  R, 


Ia4cos(40)+b4sin(40)] 


2 


[a2cos(20)+b2sin(20)l  1 


for  r  <  R2.  For  r  >  R2 


O(r,0)  =  C 


j  [a4cos(40)+b4sin(40)] 


1 


( a2cos(20)+b2sin(20)]  “ 


which  is  a  multipole  potential.  The  two  solutions  are  matched  at  R2.  We  can  now  let  R2  — > 
and  neglect  the  terms  in  Rj,  which  amount  to  a  small  correction  near  the  core,  to  get  a  potential 
valid  for  r  »  R 1 


O(r,0)  =  q 


^  [a4cos(40)+b4sin(40)l  +  ^  {a2cos(20)+b2sin(20)] 


8 


which  is  independent  of  r. 

For  numerical  calculation  we  would  need  the  piezoelectric  constant  of  the  alloy,  which 
has  not  been  determined.  To  get  an  order  of  magnitude  estimate,  we  can  use  the  piezoelectric 
constant  of  CdTe,  3  x  lO"^  C  cm*^.  The  Burgers  vector  is  4  x  10'^ cm,  and  Poisson’s  ratio  can  be 
taken  as  0.3.  Using  these  values  we  have  computed  the  charge,  potential,  and  field  plots  of 
Figure  A-2. 


13 


(a)  Piezoelectric  charge  normal  to  the  dislocation  line.  Solid  contours  are  1 . 0.5 
0.1,  and  0.05  10^^  e/cm^ going  out  radially;  dotted  contours  are  the  negative 
of  these. 


[HI] 
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(b)  Quadrupole  component  of  the  piezoelectric  potential. 


Figure  D-2.  Piezoelectric  properties  of  a  60°  dislocation  in  Hg  o.8^^o.2^^ 


(d)  Angular  variation  of  the  piezoelectric  charge  and  potential  normal 
the  dislocation  line 


Figure  D-2.  Continued 


